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Abstract: We find the high overtones of gravitational and electromagnetic quasi- 
normal spectrum of the Schwarzschild-de Sitter black hole. The calculations show 
that the real parts of the electromagnetic modes asymptotically approach zero. The 
gravitational modes show more peculiar behavior at large n: the real part oscillates 
as a function of imaginary even for very high overtones and these oscillations settles 
to some "profile" which just repeats itself with further increasing of the overtone 
number n. This lets us judge that Re(ct;) is not a constant as n ^ oo but rather 
some oscillating function. The spacing for imaginary part lm{ujn+i) — Im(ct;„) for 
electromagnetic perturbations at high n slowly approach fee as n — > oo, where kg is 
the surface gravity. In addition we find the lower QN modes for which the values 
obtained with numerical methods are in a very good agreement with those obtained 
through the 6th order WKB technique. 
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1. Introduction 

Quasinormal (QN) modes of black holes are of considerable interest recently because 
of their interpretation in ADS/CFT correspondence and of possibility to de- 

tect gravitational waves from black holes Recently it has been observed that the 
quasinormal modes can play a fundamental role in Loop Quantum Gravity (LQG) 
for asymptotically flat black holes it has been found that the asymptotic value of 
the real part of the quasinormal frequency (i.e. when the overtone approaches infln- 
ity) coincides with the Barbero-Immirizi parameter, which must be flxed to predict 
the Bekenstein-Hawking formula for entropy within the framework of LQG. All this 
stimulated development of different approaches to calculation of quasinormal modes 
r^. In particular, the analytical expression for asymptotically high QN modes of 



D-dimensional Schwarzschild black hole was obtained in [|n|] , |]T2[ . 

The asymptotic overtone quasinormal behavior was studied by NoUert |13] for 
the Schwarzschild black hole (see [5]), and in ||^ for the Schwarzschild anti-de Sitter 
black hole. Yet there is no such study for Schwarzschild-de Sitter black hole (SdS 



BH), except for near extremal case [14|. By using the 6th order WKB technique 



the low lying modes of SdS BH were estimated in |[T^. In the work |jT6|, the low 



lying QN modes for gravitational perturbations of Schwarzschild-de Sitter black hole 
were obtained with the help of the Leaver method |T^. In it was found the low 
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lying QN modes for higher (> 4) dimensional Schwarzschild black hole with different 
values of lambda term; this includes cases of Schwarzschild, Schwarzschild-de Sitter, 
and Schwarzschild-anti-de Sitter black holes. 



For the higher QN modes of the near extremal SdS BH, in |T|] it was shown that 
in contrast to asymptotically flat SBH, the QN spectrum of SdS BH have oscillatory 
behavior: the real part oscillates as a function of imaginary. Recently there appeared 
a lot of works using different numerical and analytical approaches devoted to near 
extremal Schwarzschild-de-Sitter quasinormal spectrum [pO[] . 

In this paper we make numerical investigation of the high overtones of 
Schwarzschild-de-Sitter quasinormal spectrum. We show that at very high overtones 
the real parts of the gravitational QN modes oscillates as a function of imaginary 
part. Thus the real part of u does not go to any constant limit and, thereby, the 
interpretation of the asymptotic QN frequency as those connected with the Barbero- 
Immirizi parameter in LQG is impossible. We also show that the real part of the 
electromagnetic QN modes approaches zero as an overtone number goes to infinity. 
At high overtones the QN spectrum shows rather peculiar behavior of some periodic 
weaving. The numerical results at high overtones are in a very good agreement with 
recent analytical (algebraic) equations which govern the asymptotic behavior, while 
the lower overtones are in a good agreement with those obtained earlier with the 
6-th order WKB method. 

The paper is organized as follows. In Sec. II there are basic formulas of SdS 
background for NoUert technique. Sec. Ill is devoted to the low (first ten) overtones. 
Sec. IV deals with very high overtones which let us judge about overtone asymptotic 
behavior. 



2. Basic equations 

The Schwarzschild-de-Sitter black hole is described by the metric 

ds"" = f{r)de - -7^- - r^da''- f{r) = 1 - - A^l, da^ = dO' + sin^ 9d(P', 
t[r) r 3 

(2.1) 

where M is the black hole mass, A is the cosmological constant. 

It is well known that the perturbation equations can be reduced to the 
Schrodinger wave-like equation 



d 



2 



dr*^ 

by using the so-called tortoise coordinate: 



+ uj' -V{r*)]^{r*) = 0, (2.2) 



,r- . ^. (2.3, 
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Under the choice of the positive sign of the real part of u, QNMs satisfy the 
following boundary conditions 

"^{r*) C±exp{±iujr*), r — > ±00, (2.4) 

corresponding to purely in-going waves at the event horizon and purely out-going 
waves at the cosmological horizon. 
The effective potential is given by 



V{r) = f[r) — , (2.5) 

where e = 3 for gravitational perturbations and e = for electromagnetic ones. 
The appropriate Frobenius series are 

where is the event horizon, Tc is the cosmological horizon, pe and Pc are determined 

by 

= ( - + ^— , (2.7) 



r Te / \ Tc r J \r Ve + r, 

and one can find 



Pe = 71 ^ p. 



,1 1\ /I 1 

2M + — 2M 



Substituting ( p.6|) into ( |2.2| ), we obtain the three-terms recurrent relation for a„ 



a„+ia„ + a„/?„ + cin-iln = 0, n > 0, 70 = 0, (2.8) 
where the coefficients a, (3, 7 have the form: 

rc{r, + 2re){l + n + 2pef , 2rcre(l + n + 2pe) . 
ttn = + lUJ 2.9 

n + 2pe)(n + 2pe + l)(2r2 + 2rcre-r2) ^^(rc -f Tg) 



2 ^2 



7n = J'^ , ((n + 2p,)2 - 1 - 6) 



Following Leaver [|T7| we are searching QNMs as the most stable roots of 

n Qn-l7n ^ Qn7n+1 (0 ^(W 

Pri n a„-2tn-l O «n+l7n+2 l^^.iUj 
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The infinite continued fraction on tlie riglit side of tlie equation ( 2.1CI|) converges worse 



if the imaginary part of uj increases with respect to the real part. This problem was 
circumvented by NoUert [0. He considered 

^ ^ ^2.11) 



Making use of the recurrence relation ( 2.11| ) one can find for large N: 



i?^ = Co + CiN-^'^ + C2N-^ + C^N-'^/^ . . . (2.12) 

where 



^ c e 

rJvc + 2re 



(J _ ^ / 2rgre + brcvl + 2r^ - 2(rg + rlYeiuj - A{rc + re)rlrciuJ ^q(c\ >o 



etc. 



The series (|2.12| ) converge for \uj\/N < A < oo^sowe can use this approximation 
for R]si inside the continued fraction for some N ^ — Im(u;) ~ n. In practice to find 
an appropriate we increase it until the result of the continued fraction calculations 
does not change. If one is limited by the near extremal SdS black hole the imaginary 
part, being proportional ro surface gravity, is still small in this limit, and, one can use 
the Frobenius method without NoUert modification [|T^ which includes expansion in 
A^. For non extremal values of A the situation is more complex and we have to deal 
with Nollert technique as described above. 

3. Lower overtones 

Here we present results of calculation for first ten overtones for different values of A 
and I (see Appendix in this paper). It turned out that the lower overtones obtained 
here with the help of Leaver method are in a very good agreement with those obtained 



through the 6th order WKB method Thus for example for A = 0.02 and 



/ = 2 gravitational perturbations from the 6th order WKB approach we have for the 
fundamental overtone (n = 0) = 0.3384 — 0.817i while from the continued fraction 
we obtain uo = 0.33839143 — 0.08175645i. For electromagnetic perturbations with 
/ = 1 n = we have u = 0.2259 - 0.0842i and u = 0.22594346 - 0.08410383? from 
the 6th order WKB formula and from continued fractions respectively. Note that 
the pure imaginary algebraically special value for / = 2 gravitational perturbations 
which corresponds to the 8th mode "move" to the 9th mode for A = 0.02 and to the 
higher mode for greater A. When the A is growing both the real and the imaginary 
parts of UJ are decreasing, i.e. modes damp more slowly and oscillates with greater 
real frequency. In the near extremal regime, i.e. when the A term is close to its 
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extremal value 1/9 (M = 1), the effective potential approaches the Poschl- Teller 
potential and only several first modes are well described by the formula: 



u;b = -[n + ^y+Jl{l + l)-^ (3.1) 



for scalar and electromagnetic perturbations, and by 



ub= - [n + l)i + J{l + 2){l-l)-l (3.2) 



for gravitational (axial) perturbations. Here 

54M3 



(re-3M)(rc + 6M)' 



3M. 



Note that the gravitational perturbations can be divided into the two kinds which 
can be treated separately: axial (symmetric with respect to the change —if) 



and polar. In |23| it was shown both numerically and analytically that there is the 
isospectrality between these two kinds of perturbations, i.e. they both induce the 
same QN spectrum. That is why we treat here only the axial type of gravitational 
perturbations. 



4. High overtones 

Finding of very high overtones is a time consuming procedure since the "length" 
of the continued fractions must be large enough. When we give for example the 
100000th overtone that does not mean that we found all the previous 99999 modes; 
that would require an enormous amount of time. Yet we choose those modes to 
calculate which would characterize the structure of the quasinormal spectrum at 
high overtones. 

1. Gravitational perturbations. We state that gravitational quasinormal 
spectrum of SdS black hole at asymptotically high overtones shows oscillatory be- 
havior: the real part of uo oscillates as a function of imaginary part and thus does 
not approach any constant value. We have checked this for very high overtones (see 
for example Fig. |1| where computations performed up to n ~ 165000 ). The same 
behavior for the real part of uj was observed in the near extremal regime of A . 

The imaginary part of uj is roughly proportional to n at large n and thereby can 
be approximately described by the formula: 

Im(co') ~ —ke + i, n, oo. (4.1) 
where ke is the surface gravity at the event horizon. 
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Yet, this formula is not exact even for asymptotically high overtones since the 
spacing between nearby overtones Im(ci;n+i) — Im(a;n) shows very pecuhar dependence 
on n (see as an example Fig. ^). From the first sight at figures 2 and 3 one could 
conclude that we have with numerical noise, yet if drawing the difference Im(a;„+i) — 
Im(c(j„) as a function of n in a greater scale (Fig. we see that Im(co'„_|_i) — Im(c(j„) 
shows quite ordered repeated structures thus the behavior is periodic. It is important 
that these periodic weaves do not have tendency of damping, i.e. the spacing does 
not approach ke however the average value of spacing over sufficiently large number 
of modes equals kf,: 



(Im(u;„+i) - Im(u;„)) 
N2 - iVi 



fee, n is large. (4.2) 



That is why the approximate formula ( [4.1|) is valid. The true asymptotic formula, 
apparently, must have the form: 

Im(cj) ^ + /Ce + /(n)^ , n ^ 00. (4.3) 

where f{n) is some analytically unknown part consisting of small deviations from k^. 
similar to those shown on the Fig. |^. The average value of f{n) in the sense of the 
formula ( ^^ equals zero. 

The oscillations of real part of as a function of imaginary part have very 
complicated form: the nearby overtones suffer from violate oscillations with lots of 
maximums and minimums similar to those shown on Fig. |^. These peaks form a 
larger wavy line like that shown on Fig. |l|. Thus there is little hope to find a simple 
analytical expression which could describe these oscillations. 

An important question is how we can judge whether the computed overtones 
are high enough to reflect the true asymptotic behavior?. We are sure that at a 
sufficiently large n the oscillations have a stable "profile" which just repeats itself 
with further increase of n. We can see the approaching of such a "final" profile on 
Figg. 

2. Electromagnetic perturbations. 

Electromagnetic perturbations show rather different behavior at high overtones: 
also there are oscillations of real part as a function of imaginary part, these oscilla- 
tions damp when increasing the overtone number (see Fig. ^ as an example). Thus 
the real part of u asymptotically approaches zero. We have checked already that 
at n ~ 5000 real part of uj is vanishing and we have not found QN modes with 
non- vanishing real part at higher n at all. 

The imaginary part of electromagnetic QN modes shows the same behavior at 
high overtones as the gravitational modes do. That is, even though the imaginary 
part of UJ is roughly proportional to n at large n, the spacing weaves as a function 
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0.046 - 




Figure 1: Real part of w as a function of imaginary part for some values of n = 10^, 2 • 
10"^,... In addition to this "large-scale" oscillation there is another sub-oscillations when 
considering nearby overtones (see for example Fig. 2). {1 = 2 gravitational modes, A = 
0.02) 



of n with average value Yln=Nl ^^™'^'^"+^^^"^|"^'^"^^ ke,{n is large) over sufficiently 
large quantity of modes. On contrary to gravitational perturbations, these weaves 
of the spacing of Im(u;) damp and for very high n it is approaching the equidistant 
spectrum: 

Im(co') = —ke ( n H — J i, n ^ oo. (4-4) 



5. Comparison with analytic formulas 



After the first version of this work has been appeared |25|, V.Cardoso, J.Natarrio and 
R.Schiappa ||2^ managed to find an analytical expression for n — » oo QN behavior. 
They found that the gravitational modes must obey the following algebraic equation 



Rew 




Imw 



Figure 2: Real part of u as a function of imaginary part {1 = 2 gravitational modes, 
A = 0.02). From the first sight these modes can be misunderstood as just a numerical 
noise. Yet if pictured in a wider scale like on Fig. 4 they show strict periodic structure. 

takes place for the asymptotically high (n — > oo) gravitational modes {kc is the 
surface gravity at the cosmological horizon.). While for electromagnetic modes they 
obtained the two possibilities 

uj = i{n + l/2)kf, or uj = i{n + I / 2)kc. (5.2) 

First of all the numerical data shown on Fig. ^ are in a good agreement with the 
above asymptotic formula for gravitational perturbations and the more n the closer 
numerical values to its analytical values, and, thereby, the better the QN modes 
obey the algebraic equations ^7\[ The difference is quite understood, since complete 
coincidence is only asymptotical. Yet the formula ^]T] does not predict the asymptotic 
behavior of imaginary part as a function of overtone number. It only connects the real 
and imaginary parts within one algebraic equation. From the approximate relation 

one could expect that at some very high n one could observe the equdistant 
spectrum like uj = i{n + l/2)ke. We have not observed any signs of it, since the 
weaves of the spacing of imaginary part do not show any tendency to damping. 

As to electromagnetic perturbations we observed that the asymptotic formula 
UJ = i{n + l/2)ke is really true. That is, the greater n the closer QN modes to the 
values UJ = i{n + l/2)ke- This is difficult to see for not very high overtones, but the 
coincidence of numerical and analytical results is very accurate for sufficiently high 
overtones. An instance of this approaching asymptotic regime is demonstrated on 
the following table. 
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Figure 3: Real part of a; as a function of imaginary part = 1 electromagnetic modes, 
A = 0.02). 



Comparison of numerical results shown on fig. |1| 
with the algebraic equation of ^TT] 



Numerical result 


Asymptotical formula 





039314 


- 9783.906062i 


0.037864 - 


- 9783 


907040i 





038492 


- 10228 


6244222 





037008 - 


10228 


6255192 





037073 


- 20012 


425574i 





036309 - 


20012 


4266882 





037659 


- 20457 


144537Z 





037042 - 


20457 


1455672 





038318 


- 29351 


5070862 





037857 - 


29351 


507866i 





039824 


- 30240 


943303i 





039382 - 


30240 


943943i 





039026 


- 38245 


869248Z 





038636 - 


38245 


8698692 





040054 


- 38912 


946162Z 





039657- 


38912 


9467032 





041107 


- 40024 


740368Z 





040668 - 


40024 


740824i 





041214 


- 40247 


0991582 





040764 - 


40247 


099602i 





041326 


- 40691 


816716i 





040852 - 


40691 


817141i 





041329 


- 40914 


1754912 





040843 - 


40914 


1759082 





041299 


-41136 


5342672 





040800 - 


41136 


5346762 
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Comparison of numerical results for electromagnetic / = 1 modes 
with the algebraic equation of ^]2] 



Numerical result 


Asymptotical formula 


A 


= 0.02 






analytical 




numei icai 





- 2223.5904862(n = 


10000) 


- 2223.590586i 





-4447.180971z(n = 


20000) 


-4447.181020? 





- 44471.809715i(n = 


= 200000) 


- 44471. 809720i 


A 


= 0.05 









- 1761.123825i(n = 


10000) 


- 1761.123919i 





- 3522.247649i(n = 


20000) 


- 3522.247699i 





- 35222.476494i(n = 


= 200000) 


- 35222.476499? 



At the same time we have not found any QN modes close to lu = i{n + l/2)kc 
at high n. This agrees with alternative choise of either uo = i{n + l/2)ke or uo = 
i{n + 1/2) /cc predictions. 

Thus our numerical data confirm the analytic results very well. This agreement 



between numerical and analytical results also stated in the paper 26 . 



6. Conclusion 

We have shown that even at very large overtone number the real part of the gravita- 
tional quasinormal frequency of SdS black hole oscillates as a function of imaginary 
part and these oscillations do not have any tendency to damping, i.e. the real part of 
uj asymptotically does not approach any constant value. On contrary to gravitational 
modes, real part of electromagnetic modes oscillates as well but these oscillations are 
damping with the growing of the overtone number and the real part asymptotically 
approaches zero. 

In the case of Schwarzschild black hole the interpretation of the asymptotic value 
for quasinormal frequency is known [0: the real part of it coincides with the Barbero- 
Immirizi parameter. We see that such a direct correspondence should be impossible 
for Schwarzschild de Sitter background, since the asymptotic value for quasinormal 
frequency is not a constant. Thus the possible connection of the QN frequency with 
the Barbero-Immirizi parameter in LQG is still an open question for Schwarzschild 
de Sitter black hole. 
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7. Appendix: lower overtones for electromagnetic and gravi- 
tational perturbations 



QNMs (electromagnetic) 



n 


cuM (AM2 




0.02, / = 1) 


ujM (AAf2 




0.02, / = 2) 








22594346 


— 


0.084103832 


0.41502231 


— 


0.08614382i 


1 





19922380 




0.26307549i 


0.39900686 




0.262009742 


2 





16431511 




(J.4dd81o87z 


n OTA A r\CiO A 

(J.o7(J4(Jdo4 




(J.448ozolm 


3 





13802472 




0.68534232z 


0.33705175 




0.649535542 


4 





11984831 




0.90808056z 


0.30647368 




0.862423262 


5 





10660509 




1.13166067z 


0.28135476 




1.082054132 


6 





09637285 




1.355389992 


0.26123915 




1.304842842 


7 





08827708 




1.578812402 


0.24497392 




1.528978072 


8 





08124026 




1 r\ r\c\ ^ c\ A A ^ ' 

1.802134412 


0.23156338 




1.753628192 


9 





07574049 




2.025431702 


0.22028242 




1.978417592 


10 





07054630 




2.247953382 


0.21062074 




2.203177432 


n 


uM (AM2 




0.04, / = 1) 


uM (AM^ 




0.04, / = 2) 








20061096 


— 


0.074726082 


0.36722808 


— 


0.076238762 


1 





18121373 




0.230128342 


0.35602130 




0.230646062 


2 





15204473 




0.403812122 


0.33469636 




A OA1 /(CTIIO* 

0.391451122 


3 





12839123 




0.591464162 


0.30755559 




0.562816392 


4 





11182397 




0.783988182 


0.28092911 




0.744631942 


5 





09989826 




0.976951412 


0.25829220 




0.933180152 


6 





09013754 




1.171064752 


0.23991905 




1.125126392 


7 





08368248 




1.364059792 


r\ c^cw r\r\r\ ^ o 

0.22500213 




I.3I86II3I2 


8 





07599793 




1.557181532 


0.21269507 




1.512744092 


9 





07240721 




1.751386502 


0.20235705 




1.707104442 


10 





06703704 




1.942190382 


0.19350620 




1.901499782 


n 


uM (AM^ 




0.06, I = 1) 


uM (AM2 




0.06, Z = 2) 








17089050 


— 


0.063809122 


0.31181526 


— 


0.064778182 


1 





15891839 




0.193746722 


0.30497404 




0.195150602 


2 





13694710 




0.33439045i 


0.29124711 




0.32841879i 


3 





11659604 




0.487680532 


0.27176423 




0.467631052 


4 





10180201 




0.64613954i 


0.25033424 




0.61479838i 
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Figure 4. The spacing Im(ci;„_|_i) — Im(cj„) as a function of n for gravitational perturbations, (A = 0.02, M — 1, I = 2) for large n. 
The spacing shows peculiar behavior which is NOT a numerical noise, since it is strictly periodic. 



